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Abstract 

Q-f The construction of a Monte Carlo generator for high energy hadronic and nuclear collisions 

O ■ is discussed in detail. Interactions are treated in the framework of the Reggeon Field Theory, 

D \ taking into consideration enhanced Pomeron diagrams which are resummed to all orders in the 

triple-Pomeron coupling. Soft and "semihard" contributions to the underlying parton dynamics are 
accounted for within the "semihard Pomeron" approach. The structure of cut enhanced diagrams is 
^NJ ■ analyzed; they are regrouped into a number of subclasses characterized by positively-defined contri- 

butions which dehne partial weights for various "macro-configurations" of hadronic final states. An 
iterative procedure for a Monte Carlo generation of the structure of final states is described. The 
model results for hadronic cross sections and for particle production are compared to experimental 
data. 
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1 Introduction 



Nowadays Monte Carlo (MC) generators of hadronic interactions are standard tools for data anal- 
ysis in high energy collider and cosmic ray (CR) fields. The idea behind employing such MC 
models is twofold. First of all, they provide a bridge between rigorous theoretical approaches and 
' corresponding experimental studies, thus allowing to confront novel ideas against observations. 

On the other hand, MC simulations are an inevitable part of contemporary experimental analysis 
procedures, a measurement of new phenomena depending crucially on the understanding of the 
corresponding detector response and of the contribution of the "standard" hadronic physics which 
is mimicked with the help of the MC tools. 

In particular, hadronic interaction models play an important role in investigations of very 
high energy cosmic rays. Because of the extremely low flux of such ultra-energetic particles, they 
can not be detected directly. Instead one infers their properties from measured characteristics 
of nuclear-electro- magnetic cascades, so-called extensive air showers (EAS), induced by them in 
the atmosphere. The corresponding analysis relies crucially on the MC treatment of the cascade 
development, most importantly, of its backbone - the cascade of hadron-nucleus (nucleus- nucleus) 
interactions in the atmosphere. The peculiarity of cosmic ray applications of hadronic interaction 
generators is related to the fact that one has to treat hadronic collisions at energies orders of 
magnitude higher than ones of present day colliders and that EAS characteristics depend strongly 
on model predictions for very forward spectra of secondary particles. As a consequence, CR 
interaction models, like DPMJET [1], EPOS [2], QGSJET [3], or SYBILL 0], which are designed 
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to treat general inelastic hadronic collisions, are developed in the framework of the Reggeon Field 
Theory (RFT) [5] , which allows one to take into consideration contributions from both "soft" and 
"hard" parton dynamics to the interaction mechanism. 

Soft nonperturbative interactions are described as soft Pomeron exchanges and dominate hadronic 
collisions at large impact parameters, thus giving important contributions to total, inelastic, and 
diffractive hadron-nucleus (nucleus-nucleus) cross sections. On the other hand, at sufficiently high 
energies the role of so-called semihard hadronic collisions which involve partons of moderately 
large virtualities is significantly enhanced, the smallness of the corresponding strong coupling be- 
ing compensated by large collinear and infra-red logarithms and by high density of small x partons. 
A convenient way to include such processes in the RFT treatment is provided by the "semihard 
Pomeron" approach [7] where the perturbative part of an "elementary" semihard rescattering 
is described within the DGLAP formalism, which is preceded by nonperturbative parton cascades 
("soft preevolution") described as soft Pomeron emissions. 

Additionally, high parton densities reached in "central" collisions of hadrons and, especially, 
nuclei result in significant nonlinear corrections to the interaction dynamics, related to parton 
shadowing and saturation [8]. In MC generators, such effects are typically accounted for in a phe- 
nomenological way, via energy-dependent parametrizations of some model parameters. The draw- 
back of such constructions is evident: with nonlinear effects dominating the interaction mechanism 
in the very high energy limit, model predictions are governed by the choice of the corresponding 
empirical parametrization, rather than by the underlying theoretical approach. 

In this work, we choose an alternative way, treating nonlinear interaction effects in the RFT 
framework as Pomeron- Pomeron interactions [§J [TU1 E] , based on the recent progress in the resum- 
mation of the corresponding, so-called enhanced, RFT diagrams [TSJ H31 HH ■ A MC implementation 
of such an approach has been hampered for a long time by two factors. First, with the energy 
increasing, enhanced graphs of more and more complicated topologies start to contribute signifi- 
cantly to the scattering amplitude and to partial cross sections for particular hadronic final states. 
Thus, dealing with enhanced diagrams, all-order resummation of the corresponding contributions 
is a must, both for elastic scattering diagrams and for the cut diagrams representing particular in- 
elastic processes. Secondly, it is quite nontrivial to split the complete set of cut enhanced diagrams 
into separate classes characterized by positively-defined contributions which could be interpreted 
probabilistically and employed in a MC simulation procedure. While the first problem has been 
addressed in |12l H3l [T4] , the MC implementation of the approach is discussed in the present work. 
Here we mainly address the construction of the model while the results for various particle pro- 
duction processes and applications of the model for calculations of EAS development will be the 
subject of the forthcoming publication [15] . 

The outline of the paper is as follows. In Section [21 the calculation of hadron-hadron scattering 
amplitude is discussed, taking into account enhanced diagram contributions. In Section |3l we 
consider unitarity cuts of elastic scattering diagrams and define partial contributions for various 
"macro-configurations" of the interaction, which are employed in the corresponding MC procedure, 
as described in Section |4] Finally, in Section 03 we discuss characteristic features of the developed 
model, calibration of model parameters, and present the model results for various hadronic cross 
sections. 

2 Hadron-hadron scattering amplitude 

In the RFT approach, high energy hadron-hadron scattering amplitude is defined by multiple scat- 
tering graphs of the kind depicted in Fig. [TJ The elementary rescattering contributions correspond 
to independent parton cascades developing between the projectile and target hadrons, which are 
described by Pomeron exchanges, and to parton cascades which strongly overlap in the phase space 
and interact with each other, which is described as Pomeron-Pomeron interactions. Thus, apply- 
ing the multichannel (Good- Walker-like) eikonal scheme [161 [T7] to account for contributions of 
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Figure 1: General multi-Pomeron contribution to hadron-hadron scattering amplitude; elementary 
scattering processes correspond to Pomeron exchanges (vertical thick lines) or to Pomeron-Pomeron 
interactions. 



small mass intermediate states between Pomeron emissions, elastic hadron a - hadron d scattering 
amplitude is defined as |13j 

fad(s,b)=i^ C 3/- C k/d [l-e-^W)(^)j (i) 

n ad{jk) (s, b) = 2 X v adm (s, b) + 2 X Z h w (s, b) . (2) 

Here s and b are cm. energy squared and impact parameter for the interaction, Cj/ a defines partial 
weight for hadron a elastic scattering eigenstate \j) (\a) = £\ s/Cj/a\j), J2j C 3 /a = 1)> xL(jfc) 
and X^Mjh) are ei konals corresponding to an exchange of a Pomeron or of an irreducible enhanced 
(Pomeron-Pomeron interaction) graph between the projectile and target hadrons, the latter being 
represented by eigenstates |j) and \k). 

In this work, we use the "semihard Pomeron" approach JSJ[7J to account for contributions of both 
nonperturbative soft processes and of "semihard" ones, the latter corresponding to parton cascades 
which develop at least partly in the perturbative region of relatively high virtualities |<j 2 | > Qq, 
Qq being some cutoff for pQCD being applicable. Describing the former as phenomenological soft 
Pomerons and the latter by "semihard Pomeron" exchanges, the "general Pomeron" eikonal is given 
by the sum of the two contributions: 

X V adm (s,b) = Zltm (*> b ) + X%k) b ) ■ ( 3 ) 

The soft Pomeron eikonal X^°(- f j fc ) is expressed via Pomeron emission vertices N p and the 
Pomeron propagator D r as 

Zlm^V = / d2 1 e ~ l *' b J dx l dx * Nj /a (x 1 , q 2 )Nl /d (x 2 , q 2 )D v (x 1 x 2 s,q 2 ) , (4) 

where 

D p (s,t) = 8iris (s/s ) ar e a ' lD ( s / s ^\ (5) 

with ap and a' F being the intercept and the slope of the Pomeron Regge trajectory and sq ~ 1 GeV 2 
- the hadronic mass scale. 

The Pomeron emission vertices are parametrized as 

Nj /a (x,t) = lo/a e A K« t x- a ^ (1 - x) a — , (6) 

where the exponents a par t — and ai ea d are related to intercepts of secondary Regge trajectories 

DUE]. 

The semihard contribution Xad(jk) corresponds to a piece of QCD parton ladder sandwiched 
between two soft Pomerons (see Fig. [2]) and is defined as [7J [19] 




Figure 2: A "general Pomeron" (l.h.s.) consists of the soft and semihard ones - correspondingly 
the 1st and the 2nd contributions in the r.h.s. 

Xlf m (s, b) = / A e"'* 5 J dxt dx 2 Nj /a (x u q 2 ) N r k/d (x 2 ,q 2 ) J dx+ dx~ 

x3mD r (s /x + ,q 2 )3mD p { SQ /x-,q 2 ) ^ 9i( x+ ) gj(x~) afj T> (x + x^xiX2S,Q l ) , (7) 

i,J=g,q s 

where the vertices gi for parton I (gluon or sea quark j]] coupling to the soft Pomeron are parametrized 
as 

9g{z) =r g {l-w qg ){l-zp 
g qs (z) = r g w qg I dyy^- 1 P qg (y) (1 - z/yf° , (8) 

J z 

with P qg {y) being the Altarelli-Parisi splitting function and the constant r g being fixed by momen- 
tum conservation for parton distribution functions (PDFs). 

The contribution <jj^ B (x + x~s, Qq) of parton ladder with the virtuality cutoff Qq and with the 
leg-parton types /, J and light cone momentum fractions x + , x~ is defined in a standard way (see, 

e.g. mm) 

o?? D (s,Q 2 )=Kj2 fdz + dz- [ dp* E?™(z+,QlM*)E$™,{z-,QlM*) 
i', j' J J 

x ^gpM) e(Ml-Qg), (9) 

with dajy 2 1 d-p 2 being the differential parton-parton cross section, p t - parton transverse momen- 
tum in the hard process, Afp - the factorization scale (here Mf = p\ /4), and with Ef^, (z, Qq, Q 2 ) 
describing the evolution of the parton density from the virtuality scale Qq to Q 2 . The factor 
K ~ 1.5 is designed to take effectively into account higher order QCD corrections. 

The idea behind Eq. ([7]) is to split parton evolution in an elementary scattering process in two 
parts: i) nonperturbative soft one described phenomenologically by the soft Pomeron asymptotics; 
ii) parton cascading at \q 2 \ > Qq, treated within the DGLAP formalism. The former is charac- 
terized by a significant parton diffusion in the transverse plane and, in the absence of nonlinear 
corrections, forms parton (sea quark or gluon) momentum and impact parameter distributions at 
the virtuality scale Qq HH1 |2D] . During the latter, parton transverse displacements can be ne- 
glected, leaving only the momentum-dependent part, Eq. which is characterized by a stronger 
energy-rise compared to the soft Pomeron amplitude and drives therefore the high energy behavior 
of the semihard contribution (0) . 

To calculate enhanced diagram contributions, we adopt multi-Pomeron vertices of the form |11| 

G [m ' n) = G 7 ™ +n , (10) 

1 For brevity, we do not discuss explicitly valence quark contributions to the semihard eikonal; the corresponding 
description can be found elsewhere [7l I19j . 
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Figure 3: Recursive equation for projectile net-fan contribution Xa(j)|d(fe)(y i > bi\Y, b); J/i ana - &i are 
rapidity and impact parameter distances between the projectile proton and the vertex in the handle 
of the fan. The vertex (2/2,^2) couples together 7712 projectile net-fans and n<i target net-fans. In 
addition, there are I > 1 irreducible 2-point sequences of Pomerons and Pomeron loops, exchanged 
between the vertices (yi, bi) and (y?, 62). 

where m and n are numbers of Pomerons connected to the vertex from the projectile, respectively 
target, side (m + n > 3) and the constant G is related to the triple-Pomeron coupling as 
G = r 3P /(47r 7 3). 

The eikonal X P a (j) {Vi b) f° r a "general Pomeron" exchange between hadron a (represented by 
eigenstate \ and a multi- Pomeron vertex, the two being separated from each other by rapidity 
y and transverse distance b, also receives contributions from both soft and semihard processes 

xl U) (yS = xl'$(y,b) + x F a T j) (y,b), (n) 

with the partial contributions Xa(j)* ano - ^a(j) being defined similarly to (QJ, ([7]) |20| : 

*3£M) = 87r2 ] P goe , / A e ~^" / ^1 ^ a ^i,<z 2 )^ p (^i Soe y ,<z 2 ) (12) 

= ^2 / rf 2 9 | da* Nf /a {x!,^) J dx + dx- 3mD p (s /x + ,q 2 ) 

i,J=g,q. 

where we included the vertex factors 7^ into the definition of the eikonals. 

Similarly, for a Pomeron exchange between two multi-Pomeron vertices separated by rapidity 
and impact parameter distances y and b we use 

x r (y,b) = x F °° ft (y,b) + x r ° h (y,b) (u) 

X P ^(y,b)= R J - (d\e-^ b D r {s ey, q 2 ) (15) 

07T^ I Sq e y J 

X Fsh (y,b) = ^2 / d2 1 ^ J dx+dx ~ 3mD r (s Q /x + ,q 2 ) 3mD p (s /x~, q 2 ) 

x Yl 9i(x + )gj(x-)a^ D (x + x- So ey,Q 2 ). (16) 
i,J=g,i, 

The above-defined eikonals can be used to calculate the total contribution of irreducible en- 
hanced Pomeron graphs X^m-ja- As demonstrated in [T^l the latter can be expressed via 
contributions of subgraphs of certain structure, so-called "net-fans". Those are defined by the 
Schwinger-Dyson equation of Fig. [3] and correspond to arbitrary irreducible "nets" of Pomerons, 
with neighboring net "cells" being connected to each other by 2-point sequences of Pomerons and 
Pomeron loops and with one vertex in the net having a fixed position in the rapidity and impact 
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Figure 4: Recursive representations for the contributions of irreducible 2-point sequences of 
Pomerons and Pomeron loops ^ 1oo p( 1 ) (top) and x loop (bottom), exchanged between the vertices 
(yi,6i) and (2/2,62). 



parameter spaceH The corresponding equation 



Xa«|d(fc)( 



]) (y 1 ,b 1 \Y,b) = X l a ° ( J ) (yi,b 1 ) + G / dy 2 d 2 b 2 1 - e"* ^-^l^-M) 



1-e' 



xlt ])Wk] {v2MY.b)\ e -xT ( Ua ( n( Y -y^-b- 2 \Y,b) - x ^ )|d(fe) (2/ 2 ,6 2 |r,6)l (17) 



involves the contribution X loop (2/i — 2/2, |6i — 62 1) of irreducible 2-point sequences of Pomerons and 
Pomeron loops, exchanged between the vertices (2/1,61) and (2/2,62) (2nd graph in the r.h.s. of 
Fig. [3]) , and the contribution X^alj) (2/1 > ^l) of Pomeron loop sequences exchanged between the ver- 
tex (2/1,61) and hadron a, with a single Pomeron coupled to hadron a (1st graph in the r.h.s. of 
the Figure). The 2/2 integration in the 2nd term in the r.h.s. of (|17|) is performed between £ and 
Hi — £, with £ being the minimal rapidity interval for the Pomeron asymptotics to be applicable. 
As discussed in [20], net-fan eikonals Xa(j)|d(fc) are related to parton (sea quark and gluon) distribu- 
tions which are probed during hadron-hadron interaction and are thus influenced by rescattering 
processes on the partner hadron. 

As in [14], we define the contribution x loop (2/i — 2/2, |6i — 62I) and a part of it x loop ^(yi ~~ 
2/2, 1 61 — 62 1), corresponding to Pomeron loop sequences which start from a single Pomeron coupled 
to the vertex (2/1, 61), via Schwinger-Dyson equations of Fig. 2J which gives 

X loop(1) (2/2 - 2/1, 162 - 6x1) = X P (2/2 - 2/1, |6 2 - 6i|) + G t dy' [ dH' X v (y' - y u \b' - h\) 



1-e 



-X l °°*{v2-y',\b 2 -b'\) _ x 1oop(D (j/2 _ y/} |g 2 _ g>|jl (lg) 



X loop (2/2 - 2/i, I62 - foil) = X loop(1) (2/2 - 2/i, I62 - 6i|) + G 



1-e 



-x loop (y'-yi,|fc'-bil) 



X loop (2/'-2/i,|6'-6i|) 



1-e" 



dy 1 / d 2 6' 

x loop (y2-y',|fr2-b'|) 



(19) 



In turn, for the contribution Xa°(°i > (2/i, 61) and a part of it x|°f° P ^ 1 ' > (2/i, 61), corresponding to 



2 The 2-point sequence of Pomerons and Pomeron loops, coupled to this vertex, will be referred to as the "handle 
of the fan". 
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irreducible Pomeron loop sequences with a single Pomeron coupled to the vertex (j/x, 61), this leads 
to recursive equations |14j 



x° { J)(yi,h) = x F a(j)(yi,bi) + G 



vi-£. 



dy' ld^ X : u) (y',b') 



I _ e -x loop (yi-y\\bi-b'\) _ x i°°p(i)( yi _ y \ _ b'\) 



yi-£ 



xZ (1 \yiM) = ZidyM + G J d y' j d2b ' [x:°]){y'J)-x:°t\y',b>) 

x x p (yi-y',\h-b'\). 



(20) 



(21) 



Finally, for the total contribution of irreducible enhanced graphs X^jfjk) (exchanged between 
eigenstates \j) and \k) of the projectile and of the target respectively) to elastic scattering amplitude 
one obtains 1131 [T4l 



Xlt uk) (s: b) = G dm d% { [(l - e -xS8)i*w) (l - e -*3foi»c«) - 



xTu)\ d ( k ) ~ - l/i. 16 - 61I)] [(l - e-^1.0)) e -*S8>i«« 



Xd(fe)|aO') 



+ X5 (fc) (tfi,6i) fe^ - 2/i, |6 - - Xl^\y Vx, \b~ }, (22) 

where Y = ln(s/so) and the omitted arguments of the eikonals read Xa(j)\d(k) = ^a(j)|d(fc)(^ — 
2/1,6- &x|Y,&), Xd(fe)|a(j) = Xdtk)\ a {j)(yiM\Y,b). As demonstrated in [14], Eqs. (H322D account for 
all important enhanced diagram contributions to elastic scattering amplitude. 

The generalization of the treatment for nucleus-nucleus scattering amplitude is described in 
Appendix A. 



3 Configurations of final states 



The knowledge of the elastic scattering amplitude is far insufficient for the construction of a MC 
procedure for hadronic and nuclear inelastic collisions. What we need are partial cross sections 
for particular configurations of final states, which are defined by contributions of the correspond- 
ing unitarity cuts of elastic scattering diagrams. Those can be easily derived in the nonenhanced 
eikonal scheme using the Abramovskii-Gribov-Kancheli (AGK) cutting rules [22]. Considering dia- 
grams with precisely m > 1 Pomerons being cut, each cut Pomeron corresponding to an elementary 
production process, and summing over any number of uncut Pomerons which describe additional 
elastic rescatterings, one obtains the so-called topological cross sections for hadron-hadron scat- 
tering Hang: 



( m ) / 



d 2 b ^2 Gj/aCk/d 

3,k 



(2x r ad(jk) (s,b) 



(23) 



The integrand in (|23[) can be interpreted as a probability to have precisely m elementary production 
processes in the hadron-hadron collision at impact parameter b. On the other hand, combining 
diagrams where the cut plane passes between n > 2 Pomerons, none being cut, and choosing either 
elastic or inelastic states in the cut plane for the projectile and the target, one obtains either elastic 
af d or various (low mass) diffraction cross sections. For example, for af d and for single projectile 
hadron diffraction cross section one obtains |17| : 

2 



d^b 



3,k 



1 



(24) 
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SD(proj) 



(s) = [d 2 b ^ ^ - C j/a C 1/a ) C fc/i C ra/i e ^L (3t) (^)-xL (Im) (^). (25) 



It is worth stressing that a configuration of the final state is defined by the structure of the unitarity 
cuts, here - by the number of cut Pomeron exchanges, and implies a resummation of all absorptive 
corrections due to virtual rescattering processes - uncut Pomeron exchanges. 

Taking into account enhanced Pomeron diagrams significantly complicates the analysis and 
produces a variety of final state configurations, including e.g. ones with single or multiple large 
rapidity gaps (LRG) not covered by secondary particle production. The complete set of AGK-based 
cut enhanced diagrams has been derived in [13], the corresponding contributions being composed of 
various unitarity cuts of net-fan subgraphs, and contains cut diagrams of two types. The first class 
consists of cut diagrams characterized by a "tree"-like structure of cuts; such graphs are constructed 
coupling arbitrary numbers of cut and uncut net-fan contributions in one central (not necessarily 
unique) vertex, such that each cut net-fan subgraph is characterized by a "fan'-like structure of cuts 
|13| . The diagrams of the second kind are characterized by a "zigzag"- like structure of cuts; they 
are constructed in a similar way, with the important difference that at least one of the cut net-fan 
subgraphs has a zigzaglike structure of cuts, with subsequent Pomeron end rapidities satisfying 
Vi > Vi < Vz > " • • • Treelike cut diagrams give important contributions to the total cross section 
and to partial cross sections of various final states; they provide main corrections to inclusive 
spectra of secondary particles. On the other hand, various contributions of zigzaglike cut graphs 
to the total cross section precisely cancel each other, moreover, as demonstrated in [M], they do 
not influence noticeably the rapidity gap structure of final states. Nevertheless, such diagrams 
provide contributions to inclusive particle spectra and to partial cross sections for particular final 
states. Therefore, our strategy will be to develop first a MC scheme taking into consideration 
treelike cut enhanced graphs only. After that, the procedure will be complemented by taking into 
account zigzaglike cut contributions. 

As discussed above, to obtain cross sections for various final state configurations, we shall use 
as "building blocks" contributions of various unitarity cuts of net-fan graphs. For the moment, we 
are interested in the AGK-cuts of net-fans, characterized by a fanlike structure of cuts, which are 
defined by the Schwinger-Dyson equations of Fig. [5]EI |13[ [T4] . The top line of the Figure defines 
the contribution ^iCa{j)\d(k) °^ ^ ne su bset of graphs in which the Pomeron loop sequence coupled to 
the vertex (j/i, b\) (the handle of the fan) is crossed by the cut plane. In turn, the equation in the 
bottom line gives the one of the diagrams where the handle of the fan remains uncut, 2x f a(j)\d(k) ■ 

The total contribution of fanlike cuts of net-fans is thus 2x a a (j)\d(k) = ^^a(j)\d(k) + ^^a{j)\d(k)- ^ e 
first graph in the r.h.s. of the top line corresponds to all possible AGK-based cuts of the single 
2-point sequence of Pomerons and Pomeron loops exchanged between the vertex (yi,&i) and the 
projectile hadron whereas the next two diagrams in the r.h.s. of the graphic equation describe the 
development of the cut Pomeron net: The vertex (1/2,^2) couples together m > 1 cut projectile 
net-fans, each one characterized by a fanlike structure of cuts, and any numbers m, n > of uncut 
projectile and target net-fans, such that fh+m+n > 2. There one has to subtract the contributions 
of the next two diagrams which correspond to configurations of non-AGK type, where in all the fh 
cut projectile net- fans connected to the vertex (2/27^2) the handles of the fans remain uncut and 
are situated on the same side of the cut plane, together with all the m uncut projectile net-fans. 
Finally, in the last graph in the top line of the Figure the cut plane passes between m > 2 uncut 
projectile net-fans, with at least one remaining on either side of the cut, such that a large rapidity 
gap is formed between the projectile proton and the vertex (j/2, ^2)- The diagrams in the bottom 
line of the Figure have a similar structure, as discussed in more detail in |13[ 114] , 



3 Here and in the following we use a slightly different graphic notation compared to Figs. [3] [4] a shaded ellipse with 
solid margins, positioned between the vertices (yi,6i) and (2/2, fe), denotes a general (not necessarily irreducible) 
2-point sequence of Pomerons and Pomeron loops exchanged between these vertices, with the corresponding contri- 
bution [1 — cxp(— x loop (j/i — J/2, |t>i — 62!))]- Similarly, such an ellipse with dashed margins corresponds to the AGK 
cuts of such a sequence, with the contribution 2[1 — exp(— x loop (j/i — 3/2, |6l — fe|))]- 



8 



Figure 5: Recursive representations for cut net- fan diagrams characterized by a fanlike structure 
of cuts. The top line of the figure defines the contribution 2Xa(j)\d(k) °^ ^ ne su bset of graphs in 
which the handle of the fan is cut; the bottom line gives the one of the diagrams with uncut handle, 



~fan 



2 Xa(j)\d(k)- 

As demonstrated in [13] , the total contribution of fanlike cuts of net-fans coincides with twice 
the uncut one: 

2x { :Q mk) ( yi M\Y,b) = 2xlf m{k) {yxM\YA (26) 

(j)\d(k) " 

^ )m) (yuh\Y,b) = 2 X 1 :°P(y 1 , b 1 ) + 2Gj m * ' dy 2 J <fb 2 (i _ e-x 1 ^™!^!)) 

(l - e -£7ww(y*&\ Y >^ e - 2 x3f fe )|aW) (Y-V2,b-b 2 \Y,b) _ x l : { a m(k) (y 2 , h\Y, &)] . (27) 



whereas for 2Xa(j)\d( k ) one obtains the recursive equation 



In Appendix B we derive also alternative representations for 2x^^ d ^ k y 2x^^)|d(fc) wn ich can be 
used in a MC procedure to generate the cut Pomeron structure of an irreducible cut diagram, to be 
discussed in Section [4] In addition, we obtain there subcontributions 2Xa(f)\d(k)> 2 Xa(j)\d(k)> wn ich 
correspond to such cuts of net-fan graphs in which just one cut Pomeron is coupled to hadron a. 

Using the above-defined building blocks, the complete set of cut irreducible graphs (with a 
treelike structure of cuts) for hadron- hadron scattering is given in Fig. [50 |13j . Each square bracket 
in the Figure corresponds to a positively defined contribution of a certain "macro-configuration" 
of the final state, characterized by certain topology of cut Pomerons, hence, by a definite pattern 
for secondary hadron production. For example, the diagrams in the first square bracket in Fig. [6] 
correspond to the configuration with at least two cut projectile and target net-fans (m > 2, n > 2) 
coupled together in the vertex (yi,bi), which results in a treelike structure of the final state: 
In each of the fh cut projectile net-fans any cut Pomeron may split into a few, forming a fanlike 
structure (composed of cut Pomerons) developing towards the projectile hadron (in the Os order 
with respect to the triple-Pomeron coupling a cut net-fan is represented by a single cut Pomeron 
exchanged between the hadron and the multi- Pomeron vertex); in all the h cut target net- fans 
such fanlike structures develop towards the target hadron. For the corresponding contribution we 



4 As discussed in |13| . the set of diagrams of Fig. [6] can also be represented in a form explicitly symmetric with 
respect to the projectile and the target. 

5 The 2nd and 3rd graphs in the square bracket define the subtracted contributions of the cuts of non-AGK types. 
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Figure 6: Complete set of irreducible cut diagrams characterized by a "tree"-like structure of cuts. 
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where the abbreviations are similar to the ones in 

The 2nd class of graphs corresponds to a LRG produced (in one particular inelastic rescattering 
process) between the projectile hadron and the vertex (y±, b\) and at least two cut target net- fans 
(n > 2) coupled to the vertex (j/i, bi), while in the third configuration the projectile and the target 
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interchange their places. The corresponding contributions read 
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(30) 



The next two contributions are similar to the first pair, with the difference that there is precisely 
one cut target net-fan (n = 1) coupled to the vertex (?/i,6i): 



n 
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} (31) 
•(32) 



In the 6th contribution, the secondary particles produced are separated from the target hadron 
by a LRG which extends beyond the vertex (yi,bi): In all the m > 2 cut projectile net- fans the 
handles of the fans are uncut and positioned on the same side of the cut plane, together with all 
the m > uncut projectile and n > 1 target net-fans. In the next graph the projectile and the 
target interchange their roles, the two contributions being 

rY—£ r- 

V a % k) (s,b) = 2G J d * j dH ^ \^ )ld<k) - 1 - XaWKfc)] e-xS8)i*(« (l - e-^l-o)) (33) 
rY-£ 



V a % k) (s,b)=2G I dy^tfh 



e *d(k)\a(j) — I 



-fan 



Xd{k)\a{ 3 ) 



e Xdlk)\a(j) f\ — e x£o)id(*)^. (34) 



In the graph in the 8th square bracket, there are only uncut net-fans coupled to the vertex 
(yi,6i); particle production emerges here from the cut multi-Pomeron vertex (yi,b\) only and is 
separated by large rapidity gaps from both the projectile and the target^) For the corresponding 
contribution we easily obtain 



(8) 

ad(jk) 



G 



Y-£ 



dyi Id bi 



I _ g Xa(j)|d(fc) 



I _ g Xd(fc)|a(j) 



(35) 



The next set of cut enhanced diagrams reminds the one in the 4th square bracket, being reversed 
upside-down, with the difference that the vertex (?/i, bi) is coupled to the projectile by a single cut 
sequence of Pomerons and Pomeron loops, the corresponding contributions ^Xa(^\d(k)' ^%a(j)\d(k) 
being defined in Appendix B. Thus, the cut Pomeron "tree" develops here towards the target 
while there is only one cut Pomeron coupled to the projectile. The partial contribution of such a 
configuration is 



n a % k) (s,b) = G d Vl cPh{ 



-loop -2x7u) W k) 
Xa{j)\d(k) e 



-loop e -xT( 3 )Wk) 



X a {j)\d(k) 



X 



1 - e 



gX^iaO-) _ 1 _ ^ )|a(j) 



-Xdik)\a(j) 



6 In the following we shall neglect the production of such low mass diffractive states at central rapidities, such 
that this particular set of diagrams will contribute to (quasi-)elastic rescattering processes only. As demonstrated 
in |14l . such low mass diffractive states produced at central rapidities do not provide significant contributions to 
diffraction cross sections, with the sole exception of the central diffraction (double Pomeron exchange). 
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In turn, the next contribution corresponds to a single cut sequence of Pomerons and Pomeron 
loops, exchanged between the projectile and the vertex (yi, b{) and being separated from the target 
by a LRG: 



dyi Id b\ 



-loop 

Xa(j)\d(k) e 



2 xT{i)\d{k) _ y 1o °p e -xT{i)\d{k) 



1 _ e -Xd(k)\aU) 



(37) 



Finally, the graphs in the last square bracket in Fig. [5] describe an exchange of a single cut 
sequence of Pomerons and Pomeron loops between the projectile and the target hadrons, which 
includes also a single cut Pomeron exchange, with the contribution 



O 



*a(j)\d(k) \ c 



dyi jd 2 b 1 [xd° k P )(y^ b ^ 

~ xgW) i)] + Z U) {Y - Vi, \b- 6x1) [xjfe, 61) - xS (1) (2/i, 61)] } 

As shown in [T2|, one has the identity 



(39) 



! = 1 



which relates the summary contribution of all the considered cut diagrams to the total opacity 
^lad(jk) f° r hadron-hadron scattering, Eqs. @, (|22p . Relation (|3"9")) is a direct consequence of the 
s-channel unitarity of the approach and of the fact that contributions of zigzaglike cut graphs to 
the elastic scattering amplitude precisely cancel each other [T3]. Using we can easily write 
down the absorptive cross section which corresponds to multiple secondary hadron production, 
including high mass diffraction processes: 



a*f(s)= d 2 b £ 



N 



N=l 



d 2 b 



(40) 



where the factor 



N 



ad(jk)] /N\ comes from an exchange of precisely N 



irreducible cut graphs whereas the factor exp[— Sl a d(jk)\ is obtained summing over any number 
(> 0) of elastic rescattering processes due to the exchanges of uncut graphs. Proceeding as in 
Appendix A, the treatment can be generalized to the case of nucleus-nucleus (hadron-nucleus) 
collisions, as outlined in Appendix C. 

It is noteworthy that total inelastic cross section contains also contributions from low mass 
diffraction of the projectile and/or target hadrons: 



a ad \ s ) = Vad ( s ) + Vad ( S ) + ° ad \ S ) + a ad \ s ) 



(41) 



where the latter are defined as [cf. (l2~5l) | 

^(proj)^ = J d 2 b J2 (C j/a ^-C 3/o C I/a )C fe/d C m/de 4^ t .^)-^. ( :™.M(42) 
j,k,l,m 

^ D(tai ' S) (^)= [d 2 b C i /.Q/,(%^-C i/d C iri/ ,)e-^»*)W)4V»)W)( 4 3) 



j,k,l,m 



^ad(^) = jd 2 b J2 ( C H« 6 H - Cj/aCi/a) (C k /d h m ~ C k/d C m/d ) e" («.6)-*n.-C«»0 (.,6) . (44) 

j,k,l,m 
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Eqs. (|40j|44l) form the basis for a MC treatment of inelastic hadron-hadron collisions. In par- 
ticular, using Eq. (|40p . for a given geometrical configuration of the collision (impact parame- 
ter b and elastic scattering eigenstates j and k of the projectile and target hadrons) the factor 

[O a( j(jfc)] /TV! exp[— £l a d(jk)\ can be interpreted as the probability for precisely N elementary in- 
elastic interactions to take place in the collision. Each of the elementary interactions may have 
different topologies (defined by the structure of the unitarity cuts) of the kinds discussed above, 
characterized by partial probabilities ^ad(jk) /^ad(jk) ■ It is worth stressing that the probabilistic 
interpretation of Eq. (|4"0|) and the positive-definiteness of the partial cut contributions &al(jk) 
are due to the full resummation of absorptive corrections due to virtual rescattering processes, in 
particular, due to the resummation of all the irreducible cut diagrams characterized by a given 
topology of the cuts, with any number of uncut Pomerons and any number of multi-Pomeron 
vertices. 



4 Monte Carlo procedure 

The obtained expressions allow a relatively straightforward MC implementation of the approach, 
which is discussed below for the case of hadron-hadron scattering. While low mass diffraction pro- 
cesses, being sampled according to the corresponding probabilities f??^ pr °^/ a acf 1 ' ^a^'" 8 '/ "^ 01 ' 
and cr° d D /O el [Eqs. |(I2 || -||33 || ]. are treated like in the original QGSJET model [2T] - assuming 
the PPR-asymptotics for the mass distribution of diffractive states, the "true inelastic" interactions, 
which have the partial probability cr^f/o^' , are simulated as follows. One starts from sampling 
the squared impact parameter for the collision - uniformly in the area b 2 < 6^ lax , with b max cho- 
sen sufficiently large, corresponding to negligibly small interaction probability at b > b ma x- In 
addition, one generates elastic scattering eigenstates j and k for the projectile and target hadrons 
- according to their partial weights Cj/ a) C^/d- I n the specified geometry, one defines the num- 
ber N > of elementary inelastic processes according to the Poisson distribution with the mean 
^ad(jk) ~ Eq. (|40p ; in case N = the chosen geometry is rejected and the above-discussed steps 
are repeated. 

Next, for each of the N elementary production processes one chooses first the "macro-structure" 
of the contributing cut diagrams (as defined in Fig. [5]) - according to the positively-defined weights 
ad(jk)/^ ad Q> k )'> an< ^ reconstructs the configuration of cut Pomerons for the corresponding set of 
irreducible cut graphs. For example, for the macro-configuration of the 1st square bracket in 
Fig. [5] one chooses the rapidity y\ and transverse vector b\ of the central multi-Pomeron vertex 
(yi,bi) - according to the integrand of Eq. ([2"5|) . and samples the numbers of cut projectile and 
target net-fans in, h using the Poisson distribution with the corresponding mean values 2Xafj)\d(k)' 
2Xdfk)\a(j) (rejecting the cases in, n < 2) - see the 1st term in the integrand of Eq. ([28)) . For each 
of the fh cut projectile net-fans one decides if the handle of the fan is cut - with the probability 
^a{j)\d(k) I ^a{j)\d(k)i or uncut (similarly for the h cut target net-fans); the 2nd and 3rd terms in the 
integrand of Eq. ([2"5]) are accounted for via rejection in the case all the fh cut projectile net-fans 
and/or the n cut target net- fans have their handles uncut. After that, the cut Pomeron structure 
for each of the fh + h cut net-fans is reconstructed using an iterative procedure, as discussed in 
Appendix D. By the end of the procedure one is left with cut Pomeron contributions of three types: 
i) stretched between the projectile and target hadrons; ii) between a given (projectile or target) 
hadron and a certain multi-Pomeron vertex; iii) between a pair of multi-Pomeron vertices. 

Each of those cut Pomeron contributions corresponds to an underlying elementary parton cas- 
cade developing in the respective rapidity range; hadronization of partons results in the production 
of secondary hadrons which densely fill that rapidity interval. For example, a cut Pomeron ex- 
changed between the projectile and the target hadrons gives rise to particle production in the whole 
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range [0,y]OA cut Pomeron exchanged between, say, projectile hadron and some multi-Pomeron 
vertex (y', b') results in a chain of secondaries covering the range [y 1 , Y), etc. It is noteworthy that 
we speak here about cut Pomeron contributions in the sense of Fig. [221 i-e. accounting also for 
absorptive corrections for the corresponding configuration of the final state: The cut Pomeron con- 
tribution includes also the ones of diagrams with additional multi-Pomeron vertices placed along 
the cut Pomeron line; those vertices are coupled to uncut Pomerons which are in turn connected 
to the projectile and/or the target and/or to other uncut Pomerons. 

Certain configurations obtained may contain large rapidity gaps not covered by secondary 
particles - when the corresponding rapidity intervals are not spanned by any cut Pomeron. For 
example, in the configuration of Fig. [JJ cut Pomerons cover rapidity intervals [2/1,3/2]) [j/3 ; 2/4] 5 




Figure 7 : Example of a cut enhanced graph corresponding to a final state with two LRGs (left) and 
a more complicated graph which describes absorptive corrections to the same final state (right). 

and [0,2/5], resulting in the production of chains of secondary particles in those rapidity ranges. 
Hence, LRGs are produced in the intervals [2/5,2/3] and [2/4, Y]. On the other hand, there is no 
rapidity gap in the interval [0, 2/1] which is covered by secondaries emerging from the rightmost 
cut Pomeron in the graph. While the diagram in Fig. [7J (left) is the simplest one corresponding 
to the discussed final state, the one in Fig. [JJ (right) illustrates some absorptive corrections to the 
discussed configuration, which are accounted for by the formalism. 

In addition to the already generated configuration, which is based on treelike cut enhanced 
graphs, an additional set of cut Pomerons comes from zigzaglike cut diagrams. The latter are 
treated using an effective procedure, as outlined in Appendix E. 

At the next step, both for the projectile and the target hadrons one performs energy-momentum 
sharing between all the constituent partons (Pomeron ends) connected to them and generates 
inelastic excitations of the remaining remnant states. Finally, for each cut Pomeron contribution, 
either exchanged between the projectile and the target, or between a given (projectile or target) 
hadron and a multi-Pomeron vertex, or between a pair of multi-Pomeron vertices, one chooses 
whether it is represented by its soft or semihard Pomeron component. In the latter case, like in the 
nonenhanced Pomeron scheme [3l [19] , one samples the light cone momenta for the "leg"-partons 
of the QCD ladder and performs simulation of the development of the corresponding perturbative 
parton cascade. One employs the standard treatment to reconstruct the pattern of both initial 
and final state parton emission using the forward evolution algorithms described in |19| . One 
ends up with the formation of strings stretched between the Pomeron end-point partons in case of 
soft Pomerons; for semihard Pomerons such strings are stretched also between the final s-channel 
partons resulted from the perturbative cascades, following the direction of the color flow. 

The treatment is completed with the fragmentation of strings into secondary hadrons, which is 
performed using the original procedure of the QGSJET model [21], using the algorithm described in 
|23| . with string fragmentation parameters expressed via intercepts of secondary Regge trajectories 

7 Constituent partons ("Pomeron ends") are characterized by a relatively hard light cone momentum distribution, 
hence, no LRGs arise from the energy-momentum partition between those partons and the hadron "remnant" state. 
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Figure 8: Left: Calculated total and elastic proton-proton, pion-proton, and kaon-proton cross 
sections - respectively solid, dashed, and dotted-dashed lines. Right: Calculated elastic scattering 
slope for proton-proton scattering. The compilation of experimental data (points) is from Ref. |25| . 

5 Some results and discussion 

The basic model parameters have been calibrated from the combined description of total and 
elastic hadron-proton cross sections, elastic scattering slopes, and total and diffractive structure 
functions F2, F^ 3 \ the latter two being calculated as described in [5D], generalizing the corre- 
sponding expressions to account also for Pomeron loop contributions. In turn, the parameters 
for the hadronization procedure have been tuned comparing with data on hadron production in 
proton-proton interactions, using also new data sets obtained at the Large Hadron Collider. 

Using the virtuality cutoff Qq = 3 GeV 2 between the soft and hard parton evolution, we 
obtained in particular for the soft Pomeron intercept and slope ap = 1.17, a' v = 0.11, while for 
the triple-Pomeron coupling we got r$p = 0.1 GeV, with jp = 0.4 GeV -1 . The corresponding 
results for cr^(s), OjL(s), B^,(s), d<j^ p (s,t)/dt, and for proton structure function F2(x,Q 2 ) are 
given in Figs. [8j El \W\ in comparison with experimental data. Using the simple exponential form 
(J6j) for the i-dependence of hadronic form factors, the calculated differential elastic cross sections 
agree reasonably well with measurements at small \t\ < 0.3 GeV 2 which are responsible for the 
bulk of secondary hadron production. To have a better agreement at larger values of |t|, a dipole 
parametrization for the form factor would be more suitable. 

In Fig. QTJ we compare the calculated proton diffractive structure function F^ 3 ^ (x, xp, Q 2 ) 
(proton dissociation excluded) for small Q 2 , xp, j3 = x/xp with HERA data. Our interest to this 
observable is related to its strong sensitivity to the main parameter for the enhanced Pomeron 
scheme - the triple-Pomeron coupling. It is easy to see that the model results for F^ 3 ^ agree 
with the measurements only in the limit of small Q 2 , xp, and (3. With increasing xp, the MAP 
contribution to F^ 3 ^ becomes important while for larger /3 and Q 2 so-called qq diffractive compo- 
nent (multiple Pomeron coupling to the qq-loop) has to be accounted for [35], both contributions 
neglected in the present treatment]! Thus, diffractive HERA data set the upper limit on the value 
of the triple-Pomeron vertex. 

In Fig. [T2] we plot the obtained energy dependence for single and double diffractive proton- 
proton cross sections in comparison to CDF data |36l I37| , showing also partial contributions of high 

8 Neither of the two neglected contributions involves the triple-Pomeron coupling. 
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Figure 9: Left: Calculated differential elastic proton-proton cross section for different y/s in GeV 
(as indicated in the plot). Right: Calculated differential elastic pion-proton and kaon-proton cross 
sections for the projectile lab. momentum 250 GeV/c. Experimental data are from Refs. [2"6l |2"TI 
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Figure 10: Calculated proton structure function F2(x,Q 2 ) compared to HERA data 



mass diffractive states to the discussed cross sections and the contribution of high mass diffraction of 
one proton and a low mass excitation of the other one. Here we adopt the experimental definitions 
for the diffractive cross sections, applying the respective event selection triggers to hadronic final 
states generated via a MC procedure: Single diffraction events are obtained when either a projectile 
or target proton is separated from the remaining final state of mass Mx by a LRG and M x / s < 0.15 
|36| : double diffraction events contain a central rapidity gap of size y gap > 3, which spans the 
central rapidity y — lns/2 point |37| . Diffractive states are classified as high mass ones when 
M\ > 25 GeV 2 and as low mass excitations otherwise |36| . It is noteworthy that at comparatively 
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Figure 11: Calculated proton diffractive structure function F^ 3 ^ (x, Xp, Q 2 ) compared to ZEUS 
[33] (triangles) and HI [34] (squares) data. 
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Figure 12: Calculated single <t|°(M|-/s < 0.15) and double cr^(y { g %_ 
proton cross sections (solid lines), compared to CDF data (points) \J56 



contributions to Cp°, 

proton and a low mass excitation of the other one (dotted-dashed line). 



> 3) diffraction proton- 
37] ; high mass diffraction 



<7pp (dashed lines); contribution to er^ from a high mass diffraction of one 
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low energies (y/s ~ 10 GeV) certain (theoretically) nondiffractive final states satisfy the imposed 
triggers, constituting about half of the plotted cr^p and most of the c£, D . On the other hand, 
at sufficiently high energies a part of the theoretical low mass diffraction, being described by the 
PPR-asymptotics, is classified as high mass diffraction, which explains the energy dependence of 
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the "low- high" double diffraction cross section (dotted-dashed line in Fig. H"2"T) . The obtained values 
for CTp^ agree reasonably well with the measurements, taking the fact that most of the low mass 
diffraction contribution could not be seen by the CDF detector [35]. In turn, double diffraction is 
seriously underestimated by the model. 

A comparison with selected data on secondary particle production in proton-proton and proton- 
antiproton collisions, which have been used for the model calibration, is presented in Figs.ll3 1 ll4 |ll5l 
In Figs. [13] and [14] the calculated Feynman x spectra of protons and charged pions as well as 
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Figure 13: Calculated Feynman x spectra of secondary protons in proton-proton collisions at 100 
and 200 GeV/c lab. momentum compared to data from Refs. [38] (circles) and [39] (squares). 




Figure 14: Calculated Feynman x spectra (left) and rapidity distributions (right) of positive (solid 
lines) and negative (dashed lines) pions in proton-proton collisions at 158 GeV/c lab. momentum 
compared to NA49 data [40] . 

pion rapidity distributions are plotted together with data from fixed target experiments. In turn, 
Fig. ll5l shows the results of calculations of the pseudorapidity density and of transverse momentum 
spectra of charged secondaries in non-single-diffractive (NSD) proton-antiproton and proton-proton 
collisions over a broad range of energies y/s = 0.2 7 TeV in comparison with experimental data 
from the SppS, Tevatron, and LHC colliders^ More extensive compilation of the model results for 

9 The calculations have been performed using the NSD triggers of the respective experiments. 



18 



Z 6 
-a 



- NSD 




QGSJET 11-04 


- 


✓Hr+>\i8oo 

1 JT++++T 






^200 


\ 








II 


i 1 i i i i 1 i 


i i i i i i i 


-5 


-2.5 


2.5 5 



z 



QGSJET 11-04 




1800 (|ti|< 1.0) 



> 

a 

O 



6 - 



2 - 



" NSD 




QGSJET 11-04 






KL 7000 j^<tT\ 






2360 \ 






»- > 900 -- r - f . \^ 




^* * 




iii 


i 





10 



10 



10 



-5 : 

10 r 



-5 -2.5 


2.5 5 


| NSD 
> 


QGSJET 11-04 








1 Til < 2.4 




^■*"»^ 7000 (x 10) 






= 111 


^**Vgoo (xCxrT* 



p (GeV/c) 



4 6 
p (GeV/c) 



Figure 15: Calculated pseudorapidity density (top row) and transverse momentum spectra (bottom 
row) of charged particles produced in proton-antiproton (left panels) and proton-proton (right 
panels) collisions at different cm. energies in GeV (as indicated in the plots). The data points are 
from Refs. gTJ |H S3 [H gS] . 



secondary particle production will be presented elsewhere [15] . 

An interesting potential test for hadronic interaction models, which could be performed using 
particle detectors at the LHC, has been proposed in Ref. [46]. The idea was to compare event trigger 
rates obtained by LHC experiments, using different combinations of the respective charged particle 
scintillation counters. As such counters cover a restricted range of pseudorapidities, rji < \rj\ < 772, 
with r\ x = 3.1, r? 2 = 6.5 for the TOTEM detector 07] and 771 = 2, 772 = 4 for ATLAS 03], the 
so-called minimum-bias trigger (MBT) selections by the experiments will miss a significant fraction 
of the inelastic proton-proton cross section, which will include both the low mass diffraction and a 
part of the high mass one. Using various combinations of such triggers, one gains sensitivity both 
to the absolute value of c™ el and to the contributions of single and double high mass diffraction, 
which have different selection efficiencies by such triggers. In Table [1] we present our predictions 
for trigger rates by the TOTEM and ATLAS experiments at y/s = 14 TeV, in comparison to the 
original calculations of Ref. |46| . for the selection of triggers proposed in that work: requiring 
a charge particle hit at positive rapidities only (in the interval [rji , 77 2 ] covered by scintillators) - 
MBT1, or a signal at either positive ([771,772]) or negative ([—772, —771]) rapidities - MBT2, or with 
both detectors being fired - MBT3. The calculated trigger rates for the presently attained LHC 
energy s/s = 7 TeV are collected in Table [2] Comparing our results with those of Ref. 05] , we 
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TOTEM 


ATLAS 




O"tot 


CT incl 


cmbti 


0"MBT2 


CMBT3 


cmbti 


CMBT2 


0" MBT3 


this work 


105 


76.8 


66.6 


69.7 


63.4 


63.9 


66.4 


61.3 


m 


91.5 


70.0 


50.7 


59.0 


42.4 


46.6 


50.8 


42.4 



Table 1: Calculated total, inelastic and minimum-bias (for different MBT selections) proton-proton 
cross sections at ^/s = 14 TeV. 



Ctot "mcl 


TOTEM 


ATLAS 


CMBTl CTMBT2 OMBT3 


CMBTl CTMBT2 0"MBT3 


93.3 69.7 


60.8 64.3 57.4 


58.1 60.8 55.3 



Table 2: Same as in Tabled] for V$ = 7 TeV. 



observe large differences concerning both the absolute magnitude of the calculated minimum-bias 
cross sections and for their variations between different trigger selections. Although such differences 
are partly due to hadronization effects - as we apply the respective triggers to hadronic final states 
generated via a MC procedure, the largest effect is related to a considerably higher total (hence, 
also inelastic) cross section and to a much smaller (by a factor of 2) single high mass diffraction 
cross section in our approach compared to |46| . as discussed in more detail in [TJ]. Nevertheless, 
the discussed triggers work in a similar way for both model approaches, particularly for the MBT 
selections of ATLAS: i) all the three triggers reject low mass diffraction; ii) the MBT3 trigger 
misses also most of single high mass diffraction events; iii) the MBT1 trigger rejects most of the 
target high mass diffraction. As a consequence, the ratio (ctmbt2 — ctmbt3)/(cmbt2 — ombti) 
is close to 2. Thus, the present study confirms that the method proposed in [46] will provide a 
powerful selection between different model approaches to the treatment of minimum-bias hadronic 
collisions. 

As discussed above, the model generalization to hadron-nucleus and nucleus-nucleus collisions 
does not involve any additional parameters. As an illustration, we plot in Fig. [16] the calculated 




10 

A 



Figure 16: Calculated ^4-dependence of total proton-nucleus (solid) and pion-nucleus (dashed) 
cross sections at respectively 500 and 600 GeV/c lab. momentum compared to experimental data 

m- 

A-dependence of total hadron-nucleus cross sections in comparison with experimental data. In 
addition, in Fig. [T7] the calculated Feynman x spectra and rapidity distributions of secondary 
pions in proton-carbon collisions are compared to NA49 data |50| . 

As demonstrated earlier in [12] [20], the high energy behavior of hadronic cross sections is 
strongly affected by nonlinear interaction effects. A similarly strong effect is observed when study- 



20 




Figure 17: Calculated Feynman x spectra (left) and rapidity distributions (right) of positive (solid 
lines) and negative (dashed lines) pions in proton-carbon collisions at 158 GeV/c lab. momentum 
compared to experimental data |50] , 



ing the generated configurations for hadronic final states. For proton-proton collisions, we plot in 
Fig. [18] the energy dependence of the number of "elementary pieces" of secondary production - cut 
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Figure 18: Energy dependence of the number of cut Pomerons of different types for proton- 
proton interactions: exchanged between the projectile and target protons (i), exchanged between 
the projectile or target proton and some multi-Pomeron vertex (ii), exchanged between a pair of 
multi-Pomeron vertices (iii). 

Pomerons, for the three possible contributions: i) cut Pomerons exchanged between the projectile 
and target protons - which arise from eikonal Pomeron exchanges and obtain absorptive corrections 
from various cut enhanced graphs characterized by the same pattern of the final state; ii) ones ex- 
changed between the projectile or target proton and some multi-Pomeron vertex - which come from 
cut enhanced graphs only (see the examples of the corresponding subgraphs in Fig. I23[) ; iii) cut 
Pomerons exchanged between a pair of multi-Pomeron vertices. While at moderately low energies 
the configuration of the interaction is dominated by cut Pomerons of the first type - as one would 
have in a (linear) eikonal Pomeron scheme, at higher energies the corresponding contribution is 
damped by absorptive corrections and secondary hadron production comes mainly from the other 
two contributions - which originate from cut enhanced graphs, i.e. from the treatment of nonlinear 
interactions. The physical picture behind the observed trend is obvious: The contributions of the 
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1st kind correspond to a number of elementary parton cascades developing independently between 
the projectile and target protons and hadronizing into secondary hadrons. With the energy in- 
creasing, a large number of such cascades is closely packed together in the phase space, which forces 
them to overlap and to interact with each other. Thus, production of single chains of secondaries 
covering the available rapidity range [0,Y\ is strongly reduced by virtual (elastic) rescattering of 
intermediate partons off the projectile and target hadrons. More typical become configurations of 
fan (or more complicated, see Fig. [6]) types, which correspond to multiple inelastic rescattering of 
some intermediate partons off the projectile and target hadrons, leading to a splitting or fusion 
of the cut Pomeron lines, hence, to a branching of the chains of secondaries and, generally, to a 
production of large rapidity gaps. 

Let us finally stress that the complete all-order resummation of the contributions of all sig- 
nificant enhanced diagrams performed in this work was absolutely necessary for obtaining a self- 
consistent description of hadronic cross sections and of particle production in the very high energy 
limit. This is related to the fact that in the high energy asymptotics the diagrams of the highest 
considered order with respect to Pomeron-Pomeron coupling (more precisely, the ones with the 
maximal number of Pomerons exchanged in parallel in some rapidity interval) dominate the elastic 
scattering amplitude |51| . Moreover, proceeding from one order to the next, one obtains sign- 
changing contributions. Thus, breaking the series at some given order, one obtains the total cross 
section which either falls down steeply above some energy or starts to rise in a powerlike way, vio- 
lating the unitarity bound. Hence, a meaningful answer can only be obtained after a resummation 
of the complete (infinite) series of diagrams. In fact, the situation is even more demanding when 
applied to calculations of partial cross sections for various (e.g. diffractive) hadronic final states - as 
different unitarity cuts of the same enhanced graphs give positive contributions to some processes 
while providing (negative) screening corrections to others. The simplest example of the kind is the 
triple-Pomeron diagram which provides a steeply rising contribution to high mass single diffraction 
cross section and gives rise to a strong screening correction to the single cut Pomeron cross section, 
the latter corresponding to a single elementary production process. An extensive analysis of such 
effects has been reported in the previous work |14] . where also the relative importance of various 
classes of enhanced diagrams has been investigated. 

A comment is in order on the adopted ansatz (jlOp for multi-Pomeron vertices and on the 
respective parameter 7p. As discussed already in |101 [TT] and more recently in P3], in a scheme 
based on a single Pomeron type, one is forced to choose 7p such that r^p/jp < Ap - in order 
to preserve the energy rise of the scattering amplitude. In the present scheme, apart from the 
usual soft Pomeron we have also the "semihard Pomeron" contribution [Eqs. (O, (|13p . (fl6|) ] which 
contains a perturbative "piece". In addition, we assume that multi-Pomeron coupling is dominated 
by low-g 2 processes (\q 2 \ < Qq), i.e. such vertices are coupled only to soft Pomerons or to soft "ends" 
of semihard Pomerons! 10 ! as discussed in more detail in [19l [20] . This allows us to choose 7p such 
that T3p/tp > Ap, which leads to a saturation of soft processes in the dense limit (small b and large 
s) and to a flattening of PDFs at the input scale Qq |52| . The energy rise of the scattering amplitude 
is supported at very high energies by the increase of the semihard contribution. However, neglecting 
hard (\q 2 \ > Qq) Pomeron-Pomeron coupling, we are forced to choose the Qq cutoff high enough - 
in order to safely neglect parton saturation effects at \q 2 \ > This in turn restricts our choice 

for the parameter j-p: For very small 7p, in particular, in the limit of triple-Pomeron vertices only 
(tp — > 0) , having the triple-Pomeron coupling r^p fixed by diffraction data, the saturation of the 
soft particle production would be achieved at relatively low energies over a large range of impact 
parameters. As the semihard contribution is still inefficient there, this leads to an underestimation 
of secondary particle production and to a contradiction with observations. Whether or not one can 
restrict himself with just the triple-Pomeron vertices can be investigated in the complete scheme 

10 The corresponding "hard piece" is always "sandwiched" between a pair of soft Pomerons - see Eqs. 0, 1131 1. 
1161 1 and the 2nd graph in the r.h.s. of Fig. [2] 

11 Note that for our choice of the factorization scale M| = in Eq. (|9j the chosen cutoff corresponds to the 

minimal transverse momentum in parton hard process p™ m = 2Qo — 3.4 GeV. 
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only: Taking perturbative Pomeron-Pomeron coupling into consideration. 



6 Conclusions 

In this paper, we discussed in detail the MC procedure for modeling hadronic collisions in the RFT 
framework, including the contributions of enhanced Pomeron diagrams. The principal difference 
of the presented Monte Carlo model compared to other generators of hadronic interactions is the 
direct correspondence between the RFT treatment and the MC implementation: Various hadronic 
final states are generated according to their partial cross sections. The latter are defined by the 
contributions of cut Pomeron diagrams characterized by the relevant structure of the cuts. Defining 
the contributions of certain cut subgraphs by means of recursive equations, we were able to generate 
the (generally complicated) structure of hadronic final states in an iterative fashion. 

The described model represents a self-consistent implementation of the corresponding RFT 
treatment, providing, in particular, a close link between the description of total and elastic hadron- 
proton cross sections and the generation of hadronic final states. Indeed, while the elastic scattering 
amplitude is defined by the contributions of uncut nonenhanced and enhanced diagrams, partial 
cross sections for various final states are defined by unitarity cuts of the very same diagrams, 
the summary contribution of all the cuts being related to the uncut one by s-channel unitarity 
[Eqs. (J3U); (|51))]. as demonstrated explicitly in [T3J. On the other hand, the generalization of the 
treatment to hadron-nucleus and nucleus- nucleus case, both concerning cross section calculations 
and for modeling particle production, does not involve additional adjustable parameters. 

Being based on the RFT formalism, the present treatment shares most of its usual assumptions, 
like the validity of the AGK cutting rules and eikonal vertices for Pomeron-hadron (Pomeron- 
Pomeron) coupling. It also has the usual drawback of neglecting energy-momentum correlations 
between multiple scattering processes at the amplitude level [53]. Thus, the discussed model 
remains a phenomenological one and it is experimental data which have to decide if it is suitable 
enough for the treatment of very high energy hadronic collisions. While in the current work 
we mainly addressed the construction of the MC generator and showed only some representative 
results for secondary particle spectra, a thorough comparison of the model predictions on particle 
production with available experimental data and the model applications for air shower simulation 
will be presented elsewhere |15| . 
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Appendix A 

The generalization of the approach described in Section [2] to hadron-nucleus and nucleus-nucleus 
collisions is parameter free and formally straightforward. Indeed, the only essential difference is 
that now different Pomerons in a given irreducible enhanced graph may couple to different nucleons 
of the projectile and/or target nuclei, whose positions in the impact parameter plane should be 
chosen according to the corresponding nuclear density profiles. Thus, in case of nucleus A - nucleus 
B interaction, the net-fan contribution (fl7|) should be generalized to 



m=l J i 

;-x'\%(V2,b2\Y,b)\ -xl a U(Y-y 2 ,b-b 2 \Y,b) _ not 



XA f B (yubi\Y,b)= Ex^C^I&i-^D+G^ 1 dy 2 Jd% 2 (l-e-* IO< 

XA\ B {y2M\Y,b) ,(45) 
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where y\ is the rapidity distance between the given multi-Pomeron vertex and the nucleus A, b\ 

- the position of the vertex with respect to the center of the nucleus in the transverse plane, j A 

- elastic scattering eigenstate of m-th nucleon, and b A n - its transverse vector with respect to the 
center of the nucleus. Thus, x'a\b depends on the positions {b ,b B } and the elastic scattering 
eigenstates {j A , k B } of all the A + B nucleons of the interacting nuclei, the corresponding indexes 
not shown explicitly in (|4"5")) . The general contribution of irreducible enhanced graphs is generalized 
in a similar way: 



X A n £(s,b,{j A ,k»,b A ,b B }) = G 

A 



dyi d 2 h 1 



{[(> 



not net 

Xa\b Xb a 



net 

Xa\b 



^^(Y-yiAb-h + bi 



3~Xbia I p ~Xa\ 



net 

Xb a 



E E xl ik zjyu \h - b B n \) x l ;°^{Y - yi, \b- h + bi\) - x l °°$\Y - y u \b- h + b 



(46) 



where the abbreviations are similar to the ones in 

It is obvious that Eqs. (|45fl46|) are impractical: for each particular configuration of the two 
nuclei, i.e. for each choice of the coordinates and elastic scattering eigenstates of the nucleons, one 
has to calculate Xa\b recursively, which is very time-consuming. Unlike hadron-hadron case, one 
can not make a pretabulation of the corresponding contributions which now depend on 2(A + B) 
coordinates of the nucleons, not counting the numbers of their possible eigenstates and the variables 
shown explicitly in ([43)) . 

To propose a suitable approximation for (|45fl46|) let us decompose x\ e \B as 



(47) 



with the aim to describe the dependence of X^a\ on the coordinates and eigenstates of all the 
A+B — l projectile and target nucleons but the current one, indicated symbolically by the multidot 
in the r.h.s. of (|4"T|) , by means of a single factor. Substituting ([4"7| to (|43|) and using the identity 



1 — exp 



E nct 



J2 (l-e"^)^"^ 1 ^), 



(48) 



we obtain 



m— 1 



J/l-C 



dy2 



X>^)(2/i>i&i-^i>"-) 

m— 1 

d 2 b 2 (i-e-^^-y*^-^ 



e~^k^-^ l'~>j Z$(y 2 , b 2 , Y, b, {j A , k B , b A , b B }) X ^(y 2 , \b 2 - b A \, ■ ■ ■) 

Z^(y 2 ,b 2 ,Y,b,{j A ,k B ,b A ,b B }) 

m-l B 

- E KufM & - ■ ■ ■) - E 4lk)( Y - v*> \b-h+b*\,-- •) 



1=1 



(49) 



(50) 
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Now we approximate the nuclear screening factor Z^g for m-th projectile nucleon by its value 
in the vertex (y l5 61): 

Z%l(y 2 ,b 2 , Y, b, {j A , k B ,b A , b B }) ~ Z%l( yi , h, Y, b, {j A , k B ,b A , b B }) . (51) 
Using this approximation, the solution of the nuclear net-fan equation (|45[) is 

XA f B (yi,h\Y, b) = Yl Xlf^ivi, \h - bi\, Z^l( yi M,Y, b, {j A , k B ,b A , b B })) , (52) 

m— 1 

where a partial contribution XpU) °f anv °f the A projectile nucleons is the solution of the recursive 
equation [c.f. (fT?)) ]: 



x;i} ) (yi,b',z)= x l ;^(yi,b') + G 



dy 2 J d 2 b 2 (l - e -x loop toi-w.|fc'-H) 
(1 - e-xith(y*&< z ^ Z - xJgjG,,, 6 2l Z)" 



(53) 



which can be easily tabulated as a function of its three arguments (yi, b' , Z). 
Substituting now (j5"2")) to (|4l)|) and using (|4"5)) . we obtain 



X^(s, 6, b A , b B }) =J2J2 b, {J A , k B ,b A : b B }) (54) 



m— 1 n— 1 



>Cfc, 6, {J A , k B ,b A : b B }) = G d Vl d 2 h { [( 



1 — e Xp(j 'm' I 1 1 — e * r, 'J 



E^L^XpyA) Er=l X°(fefl) nct net 



1 - e X "< fe >?> I e 



+ X P AkS) (yi, \h - b B \) [x^iY - y u \b-b 1+ bi\) - x;°^{Y - yu \b-bi + bi\)\ } , (55) 

where the omitted arguments read XpQ.^ = x£(£b)(</i, \bi~b B \, Z^ A {y u bi,Y, b, {k B ,j A , b B ,b A })), 

X P(k) = X P(k^ Y ~ yi > ^ + ^\,Z^l(Y - yi ,b-h + b A ,Y, b, {j A , k B ,b A ,b B })). 

Nucleus-nucleus elastic scattering amplitude /^s(s,6) can now be defined taking into account 
contributions from any number of Pomerons exchanged between an arbitrary pair of the projectile 
and target nucleons and from exchanges of arbitrary enhanced graphs between the two nuclei: 

fMs,b)=i((l-e-^^^S^ 



1 — exp 



I t,^AB n \sA{j A ,k B rb A ,b B }) 



m— 1 n— 1 



(56) 



^\s, b, {j A , k B ,b A , b B }) = 2 X ; pU , k s } ( Sl \b + b A ~ b B \) + 2xZnis, b, {j A , k B ,b A ,b B }) , (57) 

where for averaging over transverse coordinates and elastic scattering eigenstates of the nucleons 
we used the notation 

{h({j A , b A })) A = C itir x ■ • ■ x %/p / d%A ■ ■ ■ <*M ■■■,b A A ) h 3 A^A (b A , ...,b A ), 
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with the profile function Ta being expressed via nuclear ground state density pa as 

T A (bt ■ ■ ■ , H) = jdz? . . . dzi pAirt -.-,?£). (58) 

Expression (|57|) for nucleus-nucleus scattering amplitude reminds the usual multichannel eikonal 
form |21j . looking as a combination of binarylike nucleoli- nucleon rescatterings. In reality, each of 

the partial opacities &ab generally depends on the transverse coordinates and elastic scattering 
eigenstates of all the A + B projectile and target nucleons and contains absorptive corrections due 
to rescattering processes on those nucleons. 

Knowing the elastic amplitude, one can easily calculate total and elastic cross sections as 

*ab(«) = 2 / d 2 b3mf AB {s, b) (59) 

<b(s)= Jd 2 b\f AB (s,b)\ 2 . (60) 

Hadron-nucleus scattering amplitude is obtained using in (|55ti5T[) A = 1, 2a (&i ) = S^ib^) 
and replacing the eikonals X ¥ pp{j A kB y X^A) with X F apij A k B } , ?(£yA) for a given projectile hadron a; 
the hadron-hadron case is recovered similarly. 
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To develop a MC procedure for sampling various configurations of hadronic collisions we shall need 
alternative representations for the cut net-fan contributions 2x^ a (") \d(k)> ^Xa(j)\d(k)' Those should 
allow us to generate recursively the cut Pomeron structure for the corresponding subgraphs, with 
all the absorptive corrections due to uncut Pomerons being summed up. Such representations are 
obtained applying recursively the graphic equations of Fig. [5] to generate a i-channel sequence of 
multi- Pomeron vertices (2/2, fe), (j/3, 63), • • •, such that these vertices are coupled to uncut projectile 
and target net-fans only (mi + ri; > 1, mi = ni = 0, I = 2, 3, • • •) and they are connected to each 
other by either cut or uncut 2-point sequences of Pomerons and Pomeron loops. The recursive 
procedure stops in some vertex (y 1 , b') when: i) the vertex (y', b') is connected to the given (here, 
projectile) hadron by a cut 2-point loop sequence, with a single cut Pomeron coupled to the hadron, 
or ii) the vertex is coupled to fh! > 2 cut projectile net-fans, or hi) the vertex is coupled to m' > 2 
uncut projectile net-fans, with the cut plane positioned between them, i.e. there is a diffractive 
cut between the projectile hadron and that vertex^ The resulting Schwinger-Dyson equations are 
depicted in Figs. [T^l HOI an d read 



2xI a (- Mfc) (yi,&il^) 



2* 



loop 

a(j)\d(k) 



( yi MY,b)-2t { J )ld{k) (yiM\Y,b) 
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/ 3/1-4 
dy' I d 2 b 



E ^ ( 2 ^(i)K,)+2xi a (?)| d ( fc )) e-^i-w ~2(x^ mk) ) e-*3>.« 

n'=2 ' V 

/ \ " 

00 I Y fal \, ] 

+ (l-e-^wi'Wj -2 7 ^^ k) (y 1 ,y , M,^\Y,b)e- x ^w-^\^ ]T V ^ V 



2& mk) {y x M\Y, b) = 2x*7* (yi, h\Y, b) + G 



dy' / d 2 b 



12 Note, however, that the corresponding rapidity gap may be filled by particles resulting from other cut Pomerons 
produced in the corresponding rapidity interval. 
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Figure 19: Alternative representation for the contribution 2x^(") °f fanlike cuts of net-fans, 
the handle of the fan being cut. 




Figure 20: Alternative representation for the contribution ^XaU)\d(k) °^ f an like cuts of net-fans, 
the handle of the fan being uncut. 

x {[xjy)f3( fc ) (»!.»'. ^.^lyi^e-^i-w + x 1 °°f ) $ ik) (yi,y',b 1 ,b'\Y,b) (l-e^W^ e -^)i»( 3 ) 

oo 1 f -I -I 



m 

77i'=2 



where the arguments of the eikonals in the integrand read xT(j)\d(k) = K(j)\d(k) (f'> 6 'I F > 5 )> |d(fc) = 
X^M*)^!^. = xl a ( 5 )M(fe) (2/',?l^fc), xSw|a b ) = X%l )]a(j) (Y - y',b- b'\Y,b). 

The contributions ^Xa°j)\d(k) an( ^ ^X^(°-)|d(fc) correspond to the subset of graphs obtained in case 
(i) above, being defined by the recursive equations (see Fig. l2~Tj) 



x xZ m) {y2MYM ( e - a x»)i-«C»»Al^-x5t„^Or-»,S-&|y,8) _ x ) (63) 
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Figure 21: Recursive representations for the contributions ^Xa°f)\d(k) ~ ^Xa(f)\d(k) (^°p) ano - 
2Xa(f)\d(k) (bottom) of the subsets of fanlike cuts of net-fans, which have a single cut Pomeron 
coupled to the projectile hadron. 



ryi-Z 

xZ P mk )(yMY,b) = G dy 



, 2 d 2 b 2 (l-e-x'^iv^.r^l) 



x (l- e - x ^)W^ Y -v^-b 2 \Y,b)\ e -2 X -* )|d()e) (y 2 ,& 2 |y,&)- x '^ )|a(3) (y- a2 ,6-b 2 |yi) 
+ tm^MY^b) (e-^t^^-^.w^,i-^ _ i)]. (64 ) 

In a similar way, for the contributions Xa(j)\d(k) (V 1 ' v'> ^l; b'\Y 1 b) corresponding to t-channel 
sequences of multi-Pomeron vertices positioned between (yi,6i) and (?/,&'), coupled to uncut 
projectile and target net- fans and connected to each other and to the vertices (yi,&i), (y',b') by 
cut or uncut 2-point loop sequences [the index x (y) indicates whether the down-most (uppermost) 
loop sequence is cut, x = c (y = c), or uncut, x = u (y = u)], one obtains the equation system 



x:°^ m {yuy\bub'\YM = 

loop 



1 _ e -X lo ° P fei-?/,|bi-b'l) 



gT i dy 2 jd 2 b 2 
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I _ e -x loop (yi-a2,|fci-?2|) 
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(65) 
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y'-i 



G I dy 2 ld 2 b 2 



1 _ e -x loop (yi-2/2,|fci-&2|) 



1 _ e -X lo ° P (2/l-?A|fcl-b'l) 
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e -Xa(j)|d(fc)-^X d (fc)|a(j) + e ~ ^oWjMW X <fW|oO) — e ~ Z X a( j)|d(fc)-^X ti ( fc )| t , (i ) _J 

where the omitted arguments of the eikonals Xa(j)\d(k)> Xd(k)\a(j) are tne same as m Eqs. (|63H64p . 



Appendix C 

In case of nucleus-nucleus scattering, different cut and uncut Pomerons from the same irreducible 
graph may couple to different projectile and target nucleons. We start again from the contributions 
of fanlike cuts of net-fans, Fig. [SJ for which we obtain 
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^«5Tb(»A|v,S) _ ^ e - X sra(»A|y.?)-xSt t (y-* ! i.S-&l^ _ x^ B (y 2) 6 2 |y,6)] }, (70) 

where 2x^|^ corresponds to cut graphs where the handle of the fan is cut, 2\ l ^ B - to the ones 

where it is uncut, and ^X^b = ^^a|b + 2Xa[b i s the total contribution of fanlike cuts of net-fans; 
the dependence of the eikonals on the coordinates and eigenstates of all the nucleons is not shown 
explicitly. 

We are going to proceed like in Appendix A, expanding Xa\b^ Xa[b> ano - Xa^b as 



-fan 
XA B 



( Vl ,h\Y, b)=J2 Xfe)(tfi. |6x - bH . . .) 



m— 1 
A 



rn—1 

xZ^iyi, |6x - 6£|, . . .) = x^)(2/i, |6x - 6*|, . . .) - ^A)(yi, |6i - ft£|, . . .) 



(71) 



X%{y u h\Y, b)=J2 X^)(yi, |6x - bi\, . . .), (72) 



(73) 



and approximating the dependence on the coordinates and eigenstates of all the A+B — l projectile 
and target nucleons but the current one by some factors. 

Adding (|69|) to ([70)l . substituting the decompositions (|4"T)) . (|71ti73p . and using the identities 



e 2x pb£> 



11 e 2 ^ i4=m+lX ^V) 



(74) 
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1-exp 



m— 1 



E 



1 — e 



2 (l-e x p«4>) { 



_ y^m-l n 



J i ' e 



_ y,m-l n 



we obtain 
A 



m— 1 



m— 1 



3/1 -C 



x(l- e ^ 10 ° P ) 



2(X^,-X^ ) ) _ "2X^A ^L^te^,-^,) / ( m ) 



fan 



1 - e x pu4) ) | Z 



J A\B 



A\B 



(75) 



(76) 



(77) 



where the arguments of the eikonals in the integrand are the same as in (PI) and Z^l is defined 
in (HO]). 

Comparing ([771) to (|49|) . we obtain 



X^A)(Wi. I&1 — 6^1,...) = X^)(2/i, \h ~ bH . . .) 



(78) 



Now, substituting (jTTj) , (|71ti73|) into (JHU) , using (|74H76p , {75} , and doing the same approximation 
as in (ISTl), we obtain 



x%{yiM\Y,b) = 53 |?i - bi\,z^( yi ,b l7 Y 7 b,{j A ,k B ,b A ,b B }), 

Z ( T\yi,bi,Y, b, {j A , k B ,b A , b B }),Z B (Y - y^b-b^Y, b, {k B ,j A , b B , b A }) 
where we introduced 



Z { r\viM,Y, b, {j A , k B ,b A , b B }) = exp 



- E x^ f) (yiAh-bt\,z«> B J«\z B ) 



l=m+l 



Z A ( yi M,Y, b, {j A , k B , b A , b u }) = exp 



L 1=1 



and Xp™4-| is defined by the equation 



loop 



xfu) (yi,b',Z 1 ,Z 2 ,Z 3 ) = x l : ( °% (yi,b') + G I dy 2 \ d 2 b 2 ( 1 - e 



-X loov (vi-V2,\b'-b 2 \) 



(l - e -xlTn(y^ZuZ 2 ,Z 3 ^ Zi ^ Z3 _ -fig. ^ ^ ^ ^ ^ 



(79) 

(80) 
(81) 



(82) 



In a similar way we can obtain contributions of other subsets of cut net-fan graphs, defined in 
Appendix B, and apply them to decompose partial opacities for various macro-configurations of 
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nucleus-nucleus collisions (as defined in Fig. ^ in the form resembling binarylike nucleon-nucleon 
collisions [cf. (|M1I55|) ]: 



A B 



fiW ( a , 6, {j A , k B , b A , b B }) = E E b, {j A k B , b A , b B }) ; 



m— 1 n— 1 



which satisfy 



li 



£ n&c. 6, ^ 6 B » = n& (*> b, {j a , k B ,b A , b B }) , 



(83) 



(84) 



with being defined in (|57|) . Such a decomposition is a useful technical trick for the MC 

implementation of the approach: Each term in the decomposition corresponds to an inelastic 
rescattering process between a given pair (ran) of the projectile and target nucleons but generally 
involves additional inelastic rescatterings on other nucleons of the two nuclei. 

For example, for the configuration defined by the graphs in the 1st square bracket in Fig. [5] one 
obtains 



= d yi <Ph 



(l-e- 2 ^>) 



2x: 



2x 



n S net 



n \ — « -4 net 



P^ P {k%) _ 1 ) ^+l^(fcB) -fan 



,=m+ 



1 - e 2x p<^) ) e 



-2 E'/Ji 1 x not 



EB net 
i = lX p(fc B) 



-2 V s v nBt 
Z Xp(fcB) 13 



(85) 



Using (|84|). we can easily write down the absorptive nucleus- nucleus cross section: 



A B 

nn 



A B 



\ m,— l n—1 



E 

N mn =0 



It 11 

1 2-ft=i "«« 
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1 v mi) • 



e AB 



At B 



d % (ll- e -^=^"=^AB'^Mr\k B ,b\b B }) 



Al B 



which can be applied for a MC treatment of inelastic nucleus-nucleus (hadron-nucleus) collisions. 
It is worth stressing again that partial opacities Ornn generally correspond not only to an inelastic 
rescattering between m-th projectile and n-th target nucleons but involve also inelastic interactions 
with (m + l)-th, • • • , A-th projectile and (n + l)-th, • • • , _B-th target nucleons. 



Appendix D 

The relations of Appendix B allow one to reconstruct the cut Pomeron structure for a cut net-fan 
contribution. For simplicity, we shall illustrate the procedure neglecting the production of large 
rapidity gaps at central rapidities, i.e. neglecting the contributions of Pomeron loops and the ones 
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m'>2 



m'>2 



m'>0 



+ 




+ 




Figure 22: Recursive representation for the contribution of fanlike cuts of net-fans, Pomeron loops 
and central rapidity gaps neglected. 



-. fan 



of fanlike cuts of net-fans, which leave the handle of the fan uncut. Thus, we use ^X^j)\d(k) 
^o(i)|ci(fc)> ^a™)|d(fc) = an d the representation of Fig. [TO] takes the form (see Fig. [ 



2x t a T 3 ) Wk MMY,b) =2f a(jmk) (y 1 ,b 1 \y,b) + G / dy> / d 2 b' X ^ )|d(fc) (yi, y' , h, b'\Y, b) 



X g 2 X3(fc)| a (j) 



m'=2 



2 xl a (")|d(fc) 



2 X°U)|d(fc) _|_ ^ _ e Xa(j)|d(fc) j 



, (87) 



where the omitted arguments of the eikonals in the integrand are the same as in (|61|) and the 1st 
term in the r.h.s., ^^(j\\^rk) bi\Y, ")> ^ s * ne contribution of a i-channel sequence of cut Pomerons, 
exchanged between the vertex (yi, &i) and the projectile hadron, with the multi- Pomeron vertices 
which couple neighboring Pomerons to each other being connected to at least one uncut projectile 
or target net-fan. Similarly, 2x^ c Q\u/ k \{yi,y / ,bi,b'\Y,b) defines the contribution of such a cut 

Pomeron sequence exchanged between the vertices (2/1,61) and (y',6'). The two contributions are 
defined by recursive equations [cf. (|65[)] 

xZ(j)\d(k)(yi' b i\Y,b) = x r a ( j} (yi,bi) + G / dy 2 / d 2 b 2 x^iyi ~ 2/2, \h - b 2 \) 



x faU)\d(k)(y2,b 2 \Y,b) 



-2x:(*)| d ( fc )(y2,fc2|r,fc)-2 X ^ )|a(3) (y- !/2 ,?-h 2 |y,6) _ 1 
xlT j)]m (yi,y',h,b'\Y,b) = x r (yi - y', \h - b'\) + G f dy 2 fd 2 b 2 x r (yi - y 2 , \h - b 2 \) 



*xl? m < k) {y2,y'MAYM 



(89) 



Examples of diagrams generated by Eq. ([88)) are depicted in Fig. [23l The corresponding piece 



y,> b . 





Figure 23: Examples of cut diagrams corresponding to a single i-channel sequence of cut Pomerons 
exchanged between the vertex (yi,bi) and the projectile hadron. 



of secondary hadron production is represented by a single chain of particles produced between 
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the projectile hadron and the vertex (2/1, 61), as defined by the 1st graph in the r.h.s. of the figure 
(single cut Pomeron exchange). All the other graphs in the r.h.s. have the same particle production 
pattern and describe absorptive corrections to the process due to virtual (elastic) rescatterings on 
the projectile and target hadrons of intermediate partons of the underlying parton cascade. 

Using Eq. (|87|) . one can easily generate the cut structure for the contribution 2x ! a(j)\d(k) (fij ^ll^ 

With the probability w 1P = xZ(j)\d(k)(yi> ^\ Y > b)/x l a(j)\d(k)(yi> bi\Y, b), the particle production pat- 
tern is the one of a single cut Pomeron exchange between the projectile hadron and the vertex 
(yi, bi). In the opposite case, sampled with the probability 1 — Wxp, one generates the rapidity y' 
and transverse vector b' of the new multi-Pomeron vertex - according to the integrand of the 2nd 
term in the r.h.s. of Eq. ([87]) . Then, with the partial probability 



,-X a (j)\d(k) 



fgap 



1 - e" 



,) + ( e ^W)-l-2x^ )|d(fc) ) 



e 2 X°(j)|d(fc) 



the corresponding piece of the final state consists of a single chain of secondaries produced in 
the rapidity interval [2/1,2/] [single cut Pomeron exchange between (2/1,61) and (y',&')], with the 
projectile hadron being separated from the particles produced by a LRG. Alternatively, with the 
probability 1 — u> gap , one obtains a fanlike structure for the particle production pattern: In addition 
to the above-mentioned chain of secondaries, produced in the interval [2/1,2/'], secondary particles 
emerge from fh' > 2 cut net-fans exchanged between the vertex (2/', b') and the projectile hadron. 
In such a case, one generates the number of cut net-fans fh' according to the Poisson distribution 
with the mean ^Xa{j)\d{k) (rejecting the cases fh' < 1) and applies the above-discussed procedure 
recursively for each of the fh' cut net-fans. 

Taking into consideration Pomeron loops and central rapidity gaps, the procedure remains 
qualitatively similar, being then based on Eqs. 



(|6H68p . The difference compared to the above- 
discussed treatment is that instead of i-channel sequences of cut Pomerons (as exemplified in 
Fig. I23p one generally obtains cut i-channel sequences of Pomerons and Pomeron loops, which are 
connected to each other by multi-Pomeron vertices coupled to uncut net-fans. Hence, a similar 
algorithm is applied to reconstruct the cut Pomeron structure of those cut loop sequences. 



Appendix E 

To illustrate the effect of zigzaglike cut contributions, let us consider the simplest cut graphs of 
that kind shown in Fig. [24] (a,b). The contribution of the graph in Fig. [24] (a) is 





Figure 24: Lowest order zigzaglike cut diagrams. 



n™ d{jk) (s,b)=8G* IcfhcPh [ * dy x [ i dy 2X P aU) (Y-y 1 ,\b-bi\) 
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y-xi (l) {Y-y2Ab-b 2 \) x \ k) { yi M)x\ k) {y2M)x r {y2~ yi, \b 2 - h\), (90) 

and the one of the graph in Fig.[53](b) is defined by the same expression up to a sign, Aft^,.^ (s, b) = 
— AQ^,. fe -j(s, 6). The diagram in Fig. [24] (b) provides a (negative) screening correction to the 
eikonal configuration with two cut Pomerons. On the other hand, the one in Fig. [24] (a) introduces 
a new process, with the weight being equal to the one of the mentioned screening contribution, 
and with the particle production pattern being almost identical to the one of Fig. [24] (b); the only 
difference arises from the cut Pomeron exchanged between the vertices (yi, b\) and (jj2, b%). Thus, 
the combined effect of these two graphs is to provide additional particle production in the rapidity 
interval [yi , 2/2] ■ Hence, to account for the contributions of the graphs of Fig. [24] to secondary par- 
ticle production, one has to select final state configurations with just two cut Pomerons exchanged 
and, with the probability w zz = AfP z d(ifc) (s, b)/(ilfj^ k) (s, b) [^ 2 j^. fc) (s, 6) being the partial weight 
of the two cut Pomerons process] to add an additional cut Pomeron exchange between the vertices 
(yi,6i) and (2/2,62), with the rapidity and transverse coordinates of the vertices being generated 
according to the integrand of Eq. (|90|). At sufficiently high energies Af^^^ (s, b) > ^al(jk)( s >b) 
due to the faster energy rise of the enhanced graph contributions. A simple effective procedure 
would then be to consider w zz as the mean number of additional Pomerons to be added to the 
initial configuration. 

The general treatment of zigzaglike cut graphs follows the above-discussed logic. We re- 
strict ourselves with the set of zigzaglike cut graphs which provide nonzero contribution to in- 
clusive particle spectra and split it into two subsets whose contributions are equal up to a sign: 
K Z d(jk)( s > b ) = -^ad(jfc)( s > 6 )> where K Z d(jk) can be written as 




The particle production pattern is almost identical for the contributions fj^fjfc) an< ^ ^ad(jfe) ' exce pt 
that the former contains a cut sequence of Pomerons and Pomeron loops exchanged between the 
vertices (2/1, 61) and (2/2, 62) (with internal multi-Pomeron vertices in the sequence being generally 
coupled to uncut projectile and/or target net-fans) while the same sequence remains uncut in the 
latter contribution. Thus, the combined effect of both subsets of graphs is to add additional cut 
Pomerons resulting from this cut loop sequence. For brevity, we shall not discuss the corresponding 
technical implementation. 
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